The effect of successive summer drought periods on bacterial diversity along a plant species richness gradient

Abstract Drought is a major stressor to soil microbial communities, and the intensification of climate change is predicted to increase hydric stress worldwide in the coming decades. As a possible mitigating factor for the consequences of prolonged drought periods, above and belowground biodiversity can increase ecosystem resistance and resilience by improving metabolic redundancy and complementarity as biodiversity increases. Here, we investigated the interaction effect between plant richness and successive, simulated summer drought on soil microbial communities during a period of 9 years.To do that, we made use of a well-established biodiversity experiment (The Jena Experiment) to investigate the response of microbial richness and community composition to successive drought periods alongside a plant richness gradient, which covers 1-, 2-, 4-, 8-, 16-, and 60-species plant communities. Plots were covered from natural precipitation by installing rain shelters 6 weeks every summer. Bulk soil samples were collected 1 year after the last summer drought was simulated. Our data indicate that bacterial richness increased after successive exposure to drought, with the increase being stable along the plant richness gradient. We identified a significant effect of plant species richness on the soil microbial community composition and determined the taxa significantly impacted by drought at each plant richness level. Our data successfully demonstrates that summer drought might have a legacy effect on soil bacterial communities.


Introduction
As a result of the intensification of human activities over the past two centuries, Earth´s biosphere is facing unprecedented alterations .Climate change , resulting from increasing industrialization and air pollution, has been considered as the major driver for the change in rainfall pattern across the globe (Fowler and Hennessy 1995, Nemecek et al. 2012, Sohoulande Djebou and Singh 2016 ).As a consequence of the disturbances on the rain regime, man y r egions ar e now experiencing long periods of dr ought, suffering from insufficient precipitation or impairment of the water distribution of ri vers, lak es, and other water bodies (Eriyagama et al. 2009 ).This hydric str ess dir ectl y affects a gricultur al pr oduction, threatening the food supply chain (Osborne et al. 2013, Ostad-Ali-Askari et al. 2020, Wang et al. 2020 ), as well as the equilibrium of natural environments (Huntington 2006, Geng et al. 2015 ).
The consequences of changes in water availability can be observed on all tropic levels.Whereas data on the responses of plants and animals to w ar ds dr ought str ess ar e av ailable, still we ar e lac king a clear pictur e how micr obes in soil ar e influenced by a lack of available water, although microbes can be considered as the architects of soil quality (Vejan et al. 2016 , Docherty andGutknec ht 2019 ).To pr edict how micr obial comm unities c hange in response to drought stress is also difficult because of the complexity of the soil microbiome and the differences in ecophysiology of the single microbiota.
Bacterial comm unities r espond v ery quic kl y to envir onmental stress caused by both changes in physicochemical and biological factors, suc h as incr easing temper atur e, lac k of water, or a secondary response to a biotic factor responding to a physicochemical factor (Jansson and Hofmockel 2020 ).During water stress, bacterial cells increase intracellular solute production to achieve an osmotic equilibrium with the envir onment, whic h r equir es incr eased ener gy demands for the cell.The same happens during r e wetting, when cells r elease excess solutes into the envir onment to ac hie v e osmotic balance (Csonka 1989 ).The lack of water also reduces bacterial motility and nutrient uptake (Schimel et al. 2007 ), since the environment becomes less homogenous as water concentrations diminish.Currently, literature on the effects of drought on bacterial communities indicates that drought events are follo w ed on the one hand by an increase in Gr am-positiv e bacteria (such as Actinobacteria), which can utilize recalcitrant carbon sources and are highly present in arid, nutrient-poor soils (Connon et al. 2007).Many of these bacterial taxa are capable to generate stress-resistant structures, like spores (Zeigler 2014 ).Gr am-negativ e bacteria, on the other hand, prefer labile carbon compounds and organic nitrogen (Treseder et al. 2011), particularly in the form of plant r oot exudates, widel y abundant in eutr ophic, nutrient-ric h soils (Balasooriya et al. 2014 ).In contrast to Gr am-positiv e bacteria, r esilience to w ar ds dr ought is less pr onounced.In this sense, the intensification in both length and frequency of drought events can select drought-resistant microbes, c hanging soil micr obial population fitness and composition [15,16].T hus , the observation that drought often induces a legacy effect, e v en after r e wetting (Kaisermann et al. 2017 ), r educing soil bacterial richness in the long term is not surprising.
Taking into consideration ongoing climate change and its predicted impact on the global pr ecipitation r egime (Huntington 2006 ), investigating buffering factors of drought, is extremely relevant (Huang et al. 2023 ).One of those factors is plant community composition, which can influence and modulate the soil microbial community by recruiting and sustaining important microbial taxa (Hartman andTringe 2019 , Abedini et al. 2021 ).Pr e vious work (Wang et al. 2019, Schmid et al. 2021 ) has already shown the beneficial effects of increasing plant richness on ecosystem functions (Isbell et al. 2015 ) and microbial community composition (Lange et al. 2015, Eisenhauer et al. 2017 ).Higher plant richness promotes nutrient turnov er, biomass pr oduction, and ov er all ecosystem r esilience a gainst str ess and disturbances (Rosc her et al. 2004, Isbell et al. 2017 ).Biodiverse environments are also characterized by an increase in carbon and nitrogen stocks in soil, which ultimately contributes to higher productivity and ecosystem quality (Weisser et al. 2017, Yang et al. 2021 ).In this regard, both plant-and microbial comm unities ar e k e y to the maintenance of essential ecosystem functions , pro viding metabolic complementarity and stabilizing the ov er all ecosystem in r esponse to dr ought (Vogel et al. 2012 ).
Increasing plant richness can increase microbial activity ( Bartelt-Ryser et al. 2005, Lange et al. 2015 ) and possibly work as a mitigator of long-term drought effects on the soil microbiome.Regarding drought, plant community richness has already been shown to increase complementarity between plant species, with further adaptation to plant offspring after long periods of drought (Chen et al. 2022 ).A similar study, ho w e v er, indicates that incr easing plant diversity did not show any significant buffering effects on the soil fungal comm unity, whic h significantl y r esponded to long-term drought (Albracht et al. 2023 ).Although the effect of drought and shifts in precipitation regimes on plant communities has already been investigated (Zeppel et al. 2014 ), the complex interaction between plant diversity, microbial communities, and drought is poorly understood, especially in the phase of the recovery period after the drought event.Investigating the potential buffering effects of plant communities on soil micr obial ric hness can be crucial to maintain critical ecosystem functions as affected by climate change.
To investigate whether plant richness changes recovery of bacterial communities one year after repeated summer droughts, we made use of an experimental gradient in plant species richness, which was established in 2002 in "The Jena Experiment" (Schmid et al. 2004, Weisser et al. 2017 ).Here we simulated recurrent summer dr oughts ov er 9 y ears b y installing rain out shelters during e v ery years summer season for 6 weeks (Vogel et al. 2013a ).Control plots were also sheltered to account for potential side effects of the roof infrastructure but received ambient amounts of precipitation.Bulk soil samples were taken from control and drought treatments one year after the last summer drought treatment from all 80 plots of the plant richness gradient.Metabarcoding was used to assess changes in bacterial and arc haeal div ersity and composition across the treatments and the plant species richness gradient.Our main hypotheses were that (i) drought reduces soil micr obial ric hness and c hanges soil micr obial comm unity coposition in comparison to plots whic h r eceiv ed r egular pr ecipitation, e v en one year after the drought has been terminated and that (ii) increasing plant richness will buffer potential drought-induced soil microbial diversity loss.

Experimental design and sampling
The drought experiment was established at the field site of the Jena Experiment (Schmid et al. 2004, Weisser et al. 2017 ), which is has been running since 2002.The site is located at the floodplain of the Saale River (50 • 55 43.61 N, 11 • 35 23.64 E, altitude 130 m a.s.l.) in J ena. T he mean annual air temper atur e is 9.9 • C (1980-2010), and mean annual precipitation is 610 mm (Hoffmann et al. 2014 ).The soil is classified as an Eutric Fluvisol (World Reference Base for Soil Resources 2015 (Weisser et al. 2017 ) with a pH value range from 7.1 to 8.4 and C org 5-33 g C kg −1 (Roscher et al. 2004 ).
T he J ena Experiment consists of 80 plots (size 20 × 20 m).Those were planted with different plant community compositions in 2002.The plant communities vary in species richness (1,2,4,8,16, and 60 plant species) and number of plant functional groups (1-4 groups: grasses , small herbs , tall herbs , and legumes).Plant species richness in the manuscript is referring to this initially sown communities.Detailed information of composition and number of replicates per plant species ric hness le v el can be found in Tables S1 .The plant species richness is maintained by tw o w eeding campaigns a y ear and the plots ar e mana ged b y tw o mowing e v ents .T he plots wer e r andoml y distributed among four blocks (Fig. 1 ) to account for spatial variations in edaphic properties (including soil texture and water-holding capacity), which are related to the distance of the plots to the adjacent river Saale.For example, sand content range from 6% to 48%, silt from 38% to 71% and clay from 13% to 25% (Fischer et al. 2015 ).In this fr ame work the drought experiment was conducted from 2008 to 2016.Therefore, prior to the second annual mowing in mid-Jul y, tr anspar ent rain shelters (wood and PVC sheets) were installed for 6 weeks every year on e v ery plot (Vogel et al. 2013a ) covering two 1 × 1 m areas per plot.One roof covered the "drought" treatment, which r eceiv ed no water after installation (the 'drought' treatment), and one cov er ed the "contr ol" tr eatment, whic h r eceiv ed collected r ain water as equivalent precipitation after rain events, thereby contr olling for r oofing effects suc h as alter ed light and temper atur e (Vogel et al. 2013b ).The rain shelters reduced summer precipitation by an av er a ge of 42% in the period from 2008 to 2014 (Wagg et al. 2017 ).
As the only exception to the above described procedure, in 2013, the rain shelters were installed later (23 rd August to 23 rd September) and shorter (4, 5 weeks), because of a natural flooding e v ent of the Saale River occurring from 30 th of May to 24 th of June (Cesarz et al. 2017 ).To assess the changes in the soil micr obial comm unity after successive periods of drought, bulk soil from 0 to 15 cm was sampled in August 2017 (one year after the last experimental drought) from all plots and treatments (80 plots x 2 treatments).T herefore , 5 soil cores per treatment were pooled, resulting in approximately 50 g.The soil was sie v ed to 2 mm for homogenization and to remove bigger plant materials.Soil samples were stored at −80 • C until processing for DNA extraction and metabarcoding.

Soil moisture, pH and plant biomass
Soil moisture was measured gravimetrically.One gram of fresh soil w as w eighed (each sample w as measured in duplicate) and left in the oven overnight at 104 • C.After cooling, the soil was weighted se v er al times until a stable w eight w as r eac hed.The soil moisture is given in g of water/g of soil.pH was measured by adding 25 ml of 0.01 M calcium chloride to 10 g of air-dried soil, whic h was vigor ousl y shaken.Samples wer e incubated at r oom temper atur e for 1 h, and then the pH was measured with a calibr ated pH meter (Albr ac ht et al. 2023 ).Total biomass was calculated by harvesting the plant biomass inside the 0.1 m 2 subplot center of the 1 x 1 m plot.Samples were cut with scissors in the field at around 3 cm height above soil surface, stored in plastic bags at 4 • C for transportation to the laboratory where the total dry biomass was weighted.

DN A extr action and metabarcoding
DNA extr action fr om soil samples was performed using the DNeasy Po w ersoil Kit (Qia gen-German y).The r esulting DN A w as quantified by a Qubit fluorometric system (Thermo-Germany), and the quality was c hec ked using a Nanodrop photometric system (Thermo-Germany) and by gel electr ophor esis.As a contr ol for DN A extraction, w e included a blank extraction (DN A extraction without sample).
For metabarcoding of bacterial communities we performed Illumina next-generation amplicon sequencing targeting the 16S rRNA gene using the primer pair 515F (P ar ada et al. 2016 ) and 806R (Apprill et al. 2015 ).Each reaction had 25 μL containing 12.5 μL NEB Next High-Fidelity Master Mix (Thermo-Germany), 0.5 μL of each primer at 10 pmol/ μl, 2.5 μL of 3% BSA, 1 μl of 5 ng/ μL diluted DNA, (for the negativ e contr ol, 1 μl of DEPC-tr eated water instead) and 8 μL of DEPC-treated water.The amplification program was as follows: 98 • C for 1 min, followed by 23 cycles of 98 • C for 10 s, 55 • C for 30 s and 72 • C for 30 s, and a final extension at 72 • C for 5 min.Samples were indexed using a Nextera ® XT Index Kit v2 (Illumina-USA) and purified with Ma gSi-NGSpr ep Plus Beads (ratio 0.8 beads:1 sample) according to the manufacturer's protocol, and quality assessment was performed via a Fr a gment Anal yzer (Agilent-German y).High-quality DN A w as diluted to 4 nM and sequenced on an Illumina MiSeq using a MiSeq Reagent v3 (600 Cycle) kit.PhiX (5 pM, 20%) was loaded alongside the samples .T he ra w sequencing r eads wer e uploaded to NCBI sequencing read archive under the BioProject number PRJNA937585 and BioSample SAMN37746197.

Bioinformatics
After sequencing, samples were uploaded to the European Galaxy server ( https://usegalaxy.eu ).The Cutadpat tool was used to remov e ada pters, and r ead quality was accessed via FastQC and with the dada2 version 1.16 (Callahan et al. 2016 ) plotQualityProfile option.Trimming par ameters wer e set to 220 bp for forw ar d reads and 200 bp for reverse reads, and dada2 was used to trim the sequences without adapters.We also used dada2 to a ppl y err or r ates, mer ge the r ead pairs and make a sequence table according to the default dada2 pipeline ( https://benjjneb .github.io/dada2/tutorial.html).Taxonomy was assigned using the "assingTaxonom y and ad dSpecies" function, aligning the Amplicon Sequencing F igure 2. P earson correlation plots between alpha diversity measurements (Inverse Simpson and Observed richness) against plant species richness (natur al log).Dr ought plots ar e r epr esented in blue color, while gr ay dots and lines r epr esent the contr ol plots.Dr ought plots pr esent higher div ersity for Inverse Simpson (A), while there is no clear differentiation between the treatments in the number of observed ASVs (B).2019 ).First nonbacterial ASVs along with any ASV assigned to c hlor oplasts and mitoc hondria wer e r emov ed.Further exogenous ASVs present in the negative controls using the prevalence-based method from package decontam (v1.16) (Davis et al. 2018 ) were excluded from further analysis and batch effects from multiple sequencing runs were addressed using the "ComBat_seq" function fr om pac ka ge sv a (v3.44) (Leek et al. 2012 ).Finally, nonsingleton ASVs observed in at least 5% of the samples were used for further analysis.

Statistics and data visualization
We estimated microbial alpha diversity in each sample as richness by counting the number of ASVs and Inverse Simpson diversity index using the pac ka ges DivNet (v0.4) (Willis and Martin 2022 ) and breaka wa y (v4.7.9) (Willis and Bunge 2015 ).Differences in microbial beta diversity were estimated via PERMANOVA (999 permutations) ov er Br ay-Curtis distance matrices (Oksanen et al. 2022 ).For that, we used the adonis2 function from package Vegan v2.6.2.Additionally, we calculated the effects of the constrains soil water content, pH, plant dry biomass on community composition by using standardized estimates of beta diversity to run Redundancy analysis (RDA) using MicroViz package (v 0.10.10)(Barnett et al. 2021 ).We visualized the distance across the samples by plotting the first and second components of a singular value decomposition of the count matrix.
We analyzed the differential abundance of ASVs between control and drought samples at each plant richness level using ANCOM-BC (v1.6) (Lin and Peddada 2020 ).We filtered significantl y incr eased/decr eased ASVs ( P < 0.05), contr asting their abundances between control and drought-exposed plots across the plant richness gradient.The code and data used in these anal-yses are deposited in the GitHub repository https://github.com/rsiani/yuri _ et _ al _ 22 .
For the statistical analyses, plant species richness and number of functional groups were log-transformed to improve distribution and r eac h linearity.We fit the transformed data to a linear mixed-effects model using the lme function in the pac ka ge nlme to investigate the effects of block, plant species richness, number of functional groups, drought treatment and plot on the measured variables (Inverse Simpson and Observ ed Ric hness metrics for alpha div ersity).The dr ought tr eatment was r epr esented by subplots .T he fitting or der w as model ← lme(AlphaDiv ersityMetric ∼Bloc k + log(SpeciesRic hness) * Dr oughtTr eatment, r andom = ∼1 | Plot), moving the Block factor to a fixed term to investigate the changes in the outcome and c hanging Species Ric hness per Functional Gr oups to inv estigate the effects of functional groups instead of plant species richness.The linear mixed effect models analyzes also had the adv anta ge to account for unbalanced data (in our case the 60 species richness le v el has less observations than the other groups-Tables S1 ) by incor por ating both fixed and random effects, enabling flexible modeling of individual variability, handling unequal group sizes, and explicitly addressing the correlation structure within groups (Brown 2021 ).

Drought and plant richness effects on soil bacterial communities
Linear mixed effects model fitting r e v ealed a significant positive effect of the dr ought tr eatment and the plant ric hness gr adient on the inverse Simpson index ( P < 0.001).The number of observed ASVs (observed richness) was not significantly affected by the experimental design.The separation between drought and contr ol for inv erse Simpson can also be observ ed in the corr elation plots (Fig. 2 A), wher e dr ought-tr eated plots show higher dominance (blue dots) when compared to the control treatment (gray dots), while no clear separation could be observed in the number Table 1.(A) ANOVA table displaying the numerator and denominator degrees of freedom, alongside the F and p values for the linear mixed effects models (lme).We used Inverse Simpson and richness (number of observed ASVs) of soil bacterial communities as response variables .T he fitting order was AlphaDiv ersity ∼bloc k + log(PlantRic hness) * tr eatment, r andom = ∼1 | plot.(B) Estimated mar ginal means (EMMs) generated using the emmeans package in R. The statistical model emplo y ed w as as pr e viousl y described.Post hoc comparisons were conducted using Tuk e y's method, using pairwise comparison between control and drought across the plant diversity levels.We used the inverse Simpson and richness as microbial diversity metrics. of observed ASVs (Fig. 2 B).Ho w ever, inverse Simpson negatively correlated with plant richness ( R = −0.36 for control plots and R = −0.41 for drought plots-Fig.2 A).

Linear mixed effect
To investigate the effects of the plant richness gradient on the drought effects, we performed a post hoc test, contrasting the means of control plots against drought plots using the Inverse Simpson and number of observed ASVs as response variables.First, we observed that the difference between control and dr ought tr eated plots is significant for Inv erse Simpson along the plant richness gradient ( P < 0.001-Table 1 A), but not for richness ( P = 0.7) in alignment with the tr end observ ed in the correlation plots in Fig. 2 .Secondly, the negative values of estimated marginal means while contrasting control against drought samples presents a consistent negative result, showing that average Inv erse Simpson v alues ar e consistentl y higher in dr ought plots.We also observed that the ratio slightly decreases along the plant ric hness gr adient, suggesting that the differ ence between dr ought and control plots decreased as the plant richness increased.
To investigate the effects of the experimental design on soil bacterial composition, we used Redundancy Analysis (RDA) over Bray −Curtis dissimilarity distance (Fig. 3 ) using soil moisture, soil pH, and plant dry biomass as explanatory v ariables.We observ ed that soil bacterial community composition changes along the plant richness gradient, and PERMANOVA results indicate a significant effect of both plant species richness and drought treatments on the bacterial community composition ( P = 0.001 in both cases-Table 2 ).Although drought exerted a significant effect on soil bacterial composition, we can observe plant species richness to be the main driver of those bacterial communities (Fig. 3 ).This data also complements the r esults observ ed in the linear models and correlation plots for richness and dominance (Fig. 2 and Table 1 ), indicating that the impact of plant richness resides on the community composition, not on the overall diversity.Regarding the explanatory variables, we observed a positive correlation between the number of functional groups and plant ric hness, whic h is expected due to the increase in functional groups alongside the plant richness gradient in the Jena Experiment design.Soil moisture content displayed a slight negative correlation with the plant div ersity gr adient, while pH did not seem to be influenced by the same gradient (Fig. 3 ).

Taxonomic responses to drought and plant species richness
We also investigated the effects of both plant species richness and summer drought on the taxonomical composition of soil bacterial communities .T he taxonomical annotation of sequencing reads indicates that the overall bacterial community composition in our experiment was dominated by the same taxa, regardless of the drought treatment or plant richness level.Gemmatimonadota, Verrucomicr obiota, P atescibacteria, Myxococcota, Bacter oidota, Chlor oflexi, Acidobacteriota, Actinobacteriota, Pr oteobacteria, and Planctomycetota were the most predominant phyla, being present in all treatments and all plant richness levels (Fig. 1 ).
Differ ential abundance anal ysis (Fig. 4 ) r e v eled that monoculture plots had both increasing and decreasing ASVs belonging to the phylum Acidobacteriota with the drought treatment, both assigned to the Vicinamibacteria class, while 1 ASV assigned to the Myxococcota phylum was more abundant in the controls (Fig. 4 ).Plant richness levels 2 and 8 had both only significantly decreased ASVs after the drought treatment, assigned to Bacteroidota and Acidobacteriota phyla, respectively.On the other hand, 4-species ric hness le v el onl y pr esented a single ASV, whic h significantl y incr eased under dr ought tr eatment.This ASV was assigned to the Actinoacteriota phylum.In the 16-species plots, two ASVs wer e significantl y higher under drought and were assigned to Actinobacteriota and Bacteriodota.4 ASVs declined under drought and belonged Gemmatimonadota, Plantomycetota and Bacteroidota.The 60-species richness level was the one with F igure 3. Redundanc y analyses (RDA) plot displaying the relationships between pH, plant biomass, soil moisture, and number of plant functional groups as constrains.In the plot, the environmental variables (pH, plant biomass, soil moisture, and number of functional groups) are represented by arrows, indicating their direction and strength of influence .T he length of the arrows represents the magnitude of the effect each variable has on the biological communities .T he ellipses ar e color ed according to the div ersity le v el of eac h plot, and the dot sha pes r epr esent the dr ought tr eatment (circles for control and triangles for drought-treated plots).
Table 2. PERMANOVAs ov er Br ay −Curtis dissimilarity distance displaying the effect size and the significance of each tested variable.Plant diversity and drought treatment both presented significant effects on bacterial community composition, while the interaction between both factors was not significant.the most differential abundant ASVs with 33 ASVs.Ho w e v er, this might be an artifact caused by a lo w er number of replicated plots at this plant richness level (4 plots).The ASVs increased with the drought treatment were assigned to the Planctomycetota, Latescibacter ota, Actinobacteriota and P atescibacteria phyla, while ASVs belonging to Proteobacteriota, Planctomycetota and Chloroflexi were severely reduced.Our analyzes revealed that the affected bacterial taxa were particular to each of the analyzed plant richness levels, without a consistent bacterial taxon being affected by the drought treatment consistently.

Discussion
Our study investigated the effects of increasing plant richness on the r ecov ery of soil bacterial comm unities after r ecurr ent induced summer dr oughts.We demonstr ated that dr ought consistentl y increased bacterial dominance in comparison to plots which re-ceived ambient precipitation, even one year after the drought was terminated.This increase was consistent along the plant richness gr adient, whic h also significantl y impacted the composition of soil bacterial communities.

Bacterial di v ersity increases with long-term summer drought
The effects of drought events have already been largely studied in the soil sciences (Lipiec et al. 2013, Geng et al. 2015, Schimel 2018 ).Dr ought e v ents hav e a dr astic impact on plant and soil bacterial community abundance and composition.The lack of water incr eases temper atur e oscillations, r educes nutrient av ailability, and causes changes in the overall soil structure, making soils more compact and less porous (Lipiec et al. 2013, Zeppel et al. 2014, Geng et al. 2015, Chen et al. 2022 ).Those effects can impair plant growth and development and are largely detrimental  Manzoni et al. 2016 ), and, ob viousl y, diminishing the availability of water as a resource itself, which is essential to the basic functioning and maintenance of cellular processes (Potts 2001 ).Taking this into consideration, a drop in soil bacterial diversity shortly after a drought event is expected.Besides that, drought events can quic kl y c hange soil pr operties and nutrient av ailability, a situation wher e copiotr ophic bacterial taxa do not have enough time to react to drought exposure or simply do not possess the necessary machinery to cope with the changes imposed by drought.
In our study, ho w e v er, we tested the r ecov ery of soil bacterial community to long-term, repeated summer drought, and the effect of increasing plant species richness on this recovery.In contrast to acute drought, chronic exposure to drought stress imposes a different challenge to soil bacterial communities .T he lack of water diminishes soil homogenization, isolating bacterial communities in smaller compartments, which promotes niche formation (Carson et al. 2010 ), species differentiation (Dumbrell et al. 2010 ), and an increase in bacterial metabolic dependency (Morris et al. 2012 ).The absence of water also leaves space for more air and ther efor e mor e oxygen (Pr eece et al. 2020 ), incr easing the access of the soil bacterial communities to other gaseous and volatile substrates (Insam and Seewald 2010 ), such as methane .T he extra oxygen and new substrate availability can then be used as high energy sources for soil bacteria to explore less available and more div erse substr ates (Hartmann et al. 2011, Fest et al. 2017 ).
Our data indicate that exposure to long-term, repeated drought has a lingering effect on soil bacterial community composition.Soil samples taken 1 year after the termination of drought exposure show that diversity of soil bacterial communities increases with the drought treatment.This trend has already been observed in soils from natural holm oak forest exposed to c hr onic dr ought (Preece et al. 2019 ); ho w ever, it has-to our knowledge-never been reported in grasslands.We observed significant effects of dr ought tr eatment on the inverse Simpson (dominance) diversity index ( P < 0.001 in both cases), while richness (the number of observed ASVs) did not show any significant responses to the dr ought tr eatment ( P = 0.7).This indicates that the number of taxa did not differ between the treatments; but that the importance of dominant taxa increases in the drought treatment.This trend indicates long-lasting shifts in dominant species in each treatment, as observed in the taxonomic composition analysis.Drought is also shown to significantly affect the composition of soil bacterial comm unity, as observ ed on the PERMANOVA calculations ( P = 0.005).A possible explanation for the higher inverse Simpson index under drought might be the fr a gmentation of bacterial comm unities in the soil as a result of niche separation, alongside the promotion of less competitive bacteria in this less connected environment.Tr e v es et al. ( 2003) demonstrated that less dominant taxa have a better chance of establishing as soil moisture decreases (Tr e v es et al. 2003 ), with the competition between highly abundant taxa and less abundant taxa being more even under this condition.In this scenario, drought can reduce the nutrient availability to fast-growing taxa, allowing the growth of fastidious, less abundant taxa.Similarly, Carson et al. ( 2010 ) demonstrated that bacterial diversity increased in soils with low water content in comparison with the same soils with higher moisture (Carson et al. 2010 ).The changes in this case could be attributed to low pore connectivity due to the lack of water, increasing spatial isolation and reducing soil homogenization.In this sense, our experiment indicates a persistent difference existing in bacterial communities after 9 years of r epeated dr ought, e v en 1 year after the last tr eatment period, whic h indicates a persistent, long-term legacy effect of moder ate dr oughts on soil bacterial communities.Since plots hav e onl y been cov er ed fr om pr ecipitation during summer months, the breaks in between (other seasons when both drought exposed and control plots received equal amounts of precipitation) did not r epr esented a sufficientl y long r ecov ering time, since drought effects could still be observed.In line with our data, Vogel et al. ( 2013a ) observed a significant effect of drought treatment on the litter decomposition, also irr espectiv e of plant richness.Comparing to our context, we observed a highly significant effect of dr ought tr eatment on the soil bacterial dominance, with small correlation with the plant richness gradient.
Ho w e v er, other publications on the same experiment described contrasting effects, in comparison with our data.Albracht et al. ( 2023 ) working on the same experiment, investigated the effects of both plant ric hness incr ease and dr ought tr eatment on the diversity and composition of arbuscular mycorrhiza (AMF) and total fungal community.They reported a significant impact of the plant ric hness gr adient on the div ersity and composition of AMF and total fungal community but did not observ e an y significant effect of drought treatment on the same variables.Wagg et al. ( 2017 ) observed a reduction of plant biomass on drought treated plots in comparison to control plots, ho w ever reported a less pronounced biomass loss as the plant richness increases.Interestingl y, biomass measur ements taken in the same year of the soil sampling for this experiment (2017) did not show any significant effect of drought treatment over total plant biomass ( Fig. S2 and Table S3 ) an ymor e, indicating that r ecov ery of bacterial communities to the drought stress does not necessarily follow the patterns observed on the plant communities.Regardless of this data, we cannot completely exclude an interaction between changes in the recruitment of soil endophytic bacteria and the reduction of selectiv e pr essur e caused by the c hanges in plant-soil feedbac ks, especially since our data does not discriminate species-specific biomass or functional groups.

Plant richness gradient changes soil bacterial composition
Biodiversity is crucial for ecosystem resilience, which refers to the ability of an ecosystem to withstand and r ecov er fr om distur-bances (Cardinale et al. 2012 ).A div erse r ange of species provides functional redundancy, ensuring that ecosystem processes and services are maintained even if some species are lost (Yachi and Lor eau 1999 ).Biodiv ersity also enhances ecosystem r esistance by reducing competition through niche complementarity and incr easing ada ptability thr ough a br oader genetic pool.In this sense, plant and soil bacterial communities are intimately linked, with plant comm unities dir ectl y impacting and modulating soil bacterial communities (Liu et al. 2020 ).This modulation takes place thr ough div ersification of plant exudates in soil (Eisenhauer et al. 2017, El Moujahid et al. 2017 ), which can be used as a substrate for bacterial growth, as well as the recruitment of specific bacterial taxa to complement plant growth needs, such as phosphorus and nitrogen supply (Berihu et al. 2023 ).In drought-exposed soils, increasing plant di versity re portedly improves ecosystem resistance and resilience (Wagg et al. 2017 ), mitigating the effects of drought on biomass loss with compensatory growth after rewetting (Wagg et al. 2017 ).Ho w ever, the interaction betw een drought, the soil bacterial community and plant diversity is still poorly explored.
The effect of increasing plant richness varied according to the diversity metric analyzed.We could not observe a direct effect of plant richness on bacterial richness.We did, ho w ever, observe significant effects of plant richness gradient on the Inverse Simpson metric, as well as slightly negative correlation between this diversity metric and the plant richness gradient, which indicates that diversity calculated by Inverse Simpson slightly decreases along the plant richness gradient.This might be explained by the fact that Inverse Simpson places more emphasis on the dominance or concentration of individuals (He and Biswas 2019 ).Ther efor e, we estimate that the decrease in inverse Simpson indicated a higher proportion of dominant taxa in monocultures .T he impact of plant ric hness gr adient on the soil samples can also be observed on the PERMANOVA anal ysis, whic h shows significant effects of both plant richness gradient and drought treatment Together with the negative correlation observed between plant richness gradient, this data can indicated that increasing plant richness can favor the dominance of specific microbial species .Moreo ver, previous works on the Jena experiment observed that biodiversity effects on belowground environments might not be significant even though abov egr ound effects can be observ ed (Bessler et al. 2009 ), which can indicate that diversity effects on soil bacterial diversity can be confounded by other environmental factors.
The identity of the most abundant phyla did not change across the plant richness levels or between drought and control plots ( Fig. S1 ).Soil samples were dominated by Actinobacteria, Acidobacteria, and Proteobacteria phyla, groups commonly abundant in grassland soil samples (Janssen 2006, Fierer et al. 2007 ).The ANCOMB analyses (Fig. 4 ), ho w ever, identified ASVs which significantl y incr eased/decr eased after the dr ought tr eatment along the plant richness gradient.The selection and increase/decrease in ASV abundance does not seem to be exclusive to a single phylum, since we could observ e differ ent gener a inside of the same phylum being depleted, while others are increased.The Actinobacteriota phylum, for example, had phyla with a 2-fold increase in plant richness levels 4 and 16 (Actinomarinales order and Solirubrobacter gen us, respecti vely), while in richness level 60 two of the present ASVs were depleted ( Mycobacterium and Actinoplanes genera) while two others wer e incr eased ( Solirubrobacter genus and Micr otric hales order).Pér ez Castr o et al. ( 2019) also observ ed the decr ease in Pr oteobacteria, Verrucomicr obia, and Acidobacteria, while Actinobacteria abundance increased after dr ought str ess .T her efor e, despite the observ ation of a significant effect on the ov er all bacterial diversity, the drought treatment in our experiment did not select any specific group or taxa.
This pr efer ential accum ulation of specific taxa according to plant div ersity le v el indicates that the complementarity between bacterial and plant metabolism in the face of drought follows individual interactions at the level of bacterial species or e v en str ains l, e v en though the invers Simpson responses to drought seem to be consistent despite the plant ric hness le v els.A possible explanation for these patterns might be the variation in the number of plant functional groups implemented in parallel in the Jena Experiment (see Table S1 ) and/or the variation on species on each ric hness le v el, being the c hanges on soil bacterial composition a specie-specific interaction, instead of an overall response to the increase in surrounding diversity.The drought experiment conducted by Preece et al. ( 2019 ) also observed high variability in fungal community composition as affected by long-term drought, r epr esenting the most affected taxa highly dependent according to their experiment design.As pr e viousl y mentioned, Albr ac ht et al. ( 2023) observed the opposite: non-significant effect of drought and significant effect of plant richness gradient over fungal community.Together, this data indicates that microbial response to drought is dynamic, changes according to the investigated microbial group and according to the experimental design.

Conclusion and outlook
In summary, our data indicates that the soil bacterial community diversity was increased after long-term drought, with a rather stable response to the plant richness gradient.This response might be explained by the spatial isolation of soil bacterial communities promoted by a reduction in water potential in comparison to control samples.Changes in community composition along the plant richness gradient were observed as changes in community profile (beta diversity) instead of overall community diversity (alpha diversity) indicating that individual plant-microbe interactions might pr e v ail ov er comm unity ric hness as determining factor for soil microbial community modulation.Those findings indicate that soil bacterial diversity can adapt to long term drought conditions, being affected by the increase in plant ric hness, whic h might have important consequences for ecosystem functioning in a changing climate.In this regard, future works might address the mechanisms behind those patterns, approaching the topic with different complementary methodologies (like transcriptomic and/or metabolomic studies) and a ppr oac hing differ ent aspects of the drought stress in soil microbial communities, as the long term r ecov ery fr om str ess , for example .

Figure 1 .
Figure 1.Sc hematic r epr esentation of the Jena Experiment field site with the split-plot design regarding the drought experiment.Drought exposed and control plots are nested inside the experimental plots with increasing plant richness.Spatial distribution resembles the plot distribution in the experiment field site.

Figure 4 .
Figure 4. Differ ential anal ysis performed with ANCOMB pac ka ge .T he anal ysis estimated significant differ ences in ASV abundances between dr ought and control samples across the plant diversity gradient, with positively (right) and negatively (left) impacted taxa.Taxonomy on the left-hand side indicates the last possible le v el wher e that giv en ASV could be assigned, while bars ar e color ed according to phylum (legend in the bottom).Plant div ersity le v els ar e indicated by the numbers on the right.